Conditional Linearizability Criteria for Scalar Fourth 
Order Semi-Linear Ordinary Differential Equations 



F M Mahomed^ and A Qadir^ 

^School of Computational and Applied Mathematics, Centre for Differential Equations, 

Continuum Mechanics and Applications 

University of the Witwatersrand, Wits 2050, South Africa 

Fazal.Mahomed@wits.ac.za 

^Centre for Advanced Mathematics & Physics 

National University of Sciences & Technology 

Campus of College of Electrical Sz Mechanical Engineering 

Peshawar Road, Rawalpindi, Pakistan 

aqadirmath@yahoo . com 



Abstract. Using geometric methods for linearizing systems of second order cubically 
semi-linear ordinary differential equations and third order quintically semi-linear ordinary 
differential equations, we extend to the fourth order by differentiating the third order 
equation. This yields criteria for linearizability of a class of fourth order semi-linear 
ordinary differential equations, which have not been discussed in the literature previously. 
It is shown that the procedure can be extended to higher order. Though the results for 
the higher orders are complicated, they are doable by algebraic computing. The standard 
Lie approach, as developed at present does not seem to be amenable to giving results that 
can be handled even by algebraic computing. 

Keywords. Linearizability, symmetry algebra, geometric approach 



1 Introduction 



First order ordinary differential equations (ODEs) can always be linearized (i.e. converted 
to linear form) [1] by point transformations [2] . Lie [3] showed that all second order ODEs 
that can be converted to linear form must be cubically semi-linear 

y" + Esix,y)y" + E2ix,y)y'^ + E,{x,y)y' + Eo{x,y) = 0, (1) 

the coefficients Eq to E^ satisfying the over-determined integrable system 

12 

Ox — "g-^iy + 3"^^^ + oe — -E0-E3, 
hy = Esx - + bE2 - E.Es + eEs, 
ex = Eoy + - eEi - bEo + E0E2, 
2 1 

ey = -Ely - -E2x - be + E0E3, (2) 
1 



where b and e were (for Lie's purposes) auxiliary variables and the suffices x and y refer 
to partial derivatives. The auxilliary variables arise naturally in the geometric approach 
mentioned shortly. 

The linearization programme was carried forward to the third order by Chern [4, 5] using 
contact transformations to discuss the linearizability of equations reducible to the linear 
forms u"'{t) — and u"'{t) + u{t) — and Grebot [6, 7] used point transformations that 
were restricted to the class of transformations t = 0(x), u = ilj{t, x) for the same purpose. 
It was subsequently shown [8] that there are three classes of third order ODES that are 
linearizable by point transformations, viz. those that reduce to the above two forms or 
u"'{t) + a{t)u{t) = 0. (Note that one cannot choose a to be a constant, as it would 
then reduce to the second form by re-scaling the independent variable.) After that Neut 
and Petitot [9] and later Ibragimov and Maleshko [10] used the original Lie procedure 
[3], of postulating some point transformation that yields the desired linearizability, and 
determining the criteria for it to exist. The resulting criteria are very complicated and 
barely manageable for algebraic computing. Clearly, higher order attempts would not be 
feasible for this method. Maleshko [11] also provided a simple algorithm to reduce third 
order odes of the form y'" = F[y,y',y"] to second order and used the d/dt symmetry 
to reduce the order to two. He then used the Lie linearizability criteria to determine 
the linearizability of the third order equation. This could be trivially extended to the 
fourth order for y"" = F[y' ,y" ,y"'], but this extension would not work for any non-trivial 
problem. 

Though no progress was made towards providing explicit procedures for linearizing sys- 
tems of ODEs, it was shown [12] that there will generally be multiple classes, even for 
second order systems. The number of classes for even dimensions, 2m, is 2m^ + 3 and for 
odd dimensions, (2m — 1) it is 2m^ — 2m -|- 4. Clearly, there was no question of proceeding 
to higher order systems with the previous procedures. 

Using the connection between the isomctry algebra and the symmetries of the system of 
geodesic equations [13, 14], UncarizabiUty criteria were stated for a system of second order 
quadratically semi-hnear ODEs, of a class that could be regarded as a system of geodesic 
equations [15], that we will call of geodesic type. The criteria came from requiring that 
the coefficients in the equations, regarded as Christoffel symbols, yield a zero curvature 
tensor. The criteria also provide the required compatibility conditions. Further, even 
for larger dimensional systems, it is possible to associate a metric with the system of 
geodesies when these criteria are met [16]. The flatness of the metric allows a coordinate 
transformation to be deflned from the given metric to a Cartesian (Eulidean or pseudo- 
Euclidean) form. These are the linearizing transformations. They can be determined 
using complex variables [15] and used to write down the solution in terms of the original 
variables. 

Utilizing the projection procedure of Aminova and Aminov [13], which appeals to the fact 

that the geodesic equations do not depend on the geodetic parameter, the system of n 
second order ODEs of geodesic type can be reduced to system of (n — 1) second order 



cubically semi-linear ODEs [17]. When applied to a system of two dimensions we get a 
scalar cubically semi-linear ode that yields the Lie cirteria! Apphed to a system of three 
dimensions, one obtains a system of two cubically semi-linear odes with extended Lie 
cirteria. The system of odes so obtained is in the class of maximally symmetric equations 
(out of the total of five classes mentioned earlier) . 

Differentiating the quadratically and cubically semi-hnear system of ODEs relative to the 
independent variable gives third order ODEs. Taking the general class of the scalar third 
order ODE one gets linearizability criteria for scalar third order ODEs [18]. This class is is 
not included in the Neut and Petitot [9] and Ibragimov and Maleshko classes [10]. Though 
there can be an overlap with the Maleshko class [11] it is obviously not contained in that 
either. It may be wondered how it can be outside the three classes allowed by [8]. The 
reason is that it is not obtained by point transformations but is linearized conditionally to 
a second order ODE being linearizable. In effect its first integral is a linearizable second 
order semi-linear ODE. This is, thus, a new type of linearizability. 

The quadratically semi-linear system of geodesic type, on differentiation yields a cubically 
semi-linear third order system of the same dimension. This is obtained by replacing the 
second order term in the differentiated equation by a quadratic expression using the 
original equation. It could also be written as a third order ODE with a second order 
term multiplied by a first order term and then followed by a term quadratic in the first 
derivative. The linearizability criteria for this system for two dimensions have also been 
derived [19]. 

It is clear that the geometrical approach has far-reaching consequences for the lineariz- 
ability of ODEs. In this paper it will be used for stating criteria for the conditional 
linearizability of fourth order ODEs in some detail aand its utilization for higher orders 
will also be discussed. First a brief summary of the notation used, relevant for our present 
purposes, will be provided. 



2 Notation and review 

Equations of geodesic type for two dimensions are of the form 



There are six coefficients and they may be identified with the six Christoffel symbols for 
a 2-dimensional space. 



x" = a(a:, y)x''^ + 2b{x, y)x'y' + c(a;, y)y''^i 
y" = d{x, y)x''^ + 2e{x, y)x'y' + f{x, y)y''^. 



(3) 




(4) 



Q 



Requiring that the space corresponding to these Christoffel symbols is flat 

pi pi pi I pi pm pi pm a /c \ 

^jkl ~ ^ jl,k ^ jk,l ^ mfe^ ^ mi^ jfe " W) 

yields the linearizability conditions 

tty — bx + be — cd = 0, by — + [ac — b^) + {bf — ce) = 0, 
dy-e^- {ae - bd) - {df - e^) = 0, (6 + f)^ = (a + e)y, (6) 

which provide the metric coefficients through 

Px — —2{ap + dq), = —bp — (a + e)q — dr, = —2{bq + er), 
Py = -2{bp + eq), qy = -cp - (6 + f)q - er, Ty = -2{cq + fr). (7) 

The compatibility of these equations is guaranteed by the curvature tensor being zero. 
The metric coefficients being obtained the linearizing transformation is available as a 
coordinate transformation and its inverse can then be computed. 

The above system of two equations can be projected to the scalar cubically semi-linear 
ODE 

y" + cy" -if- 2b)y" + {a- 2e)y' - d = 0, (8) 

which is clearly of the form (1) for suitable values of the coefficients in the former equation. 

It is here that one sees the "auxilliary variables" arising naturally. Note that there arc only 
four coefficients here as against the six Christoffel symbols and for the original system of 
two equations. This degeneracy allows for extra classes when this procedure is applied to 
larger systems of ODEs but does not interfere with the uniqueness for the scalar equation 
[16]. We shall write g — f — 2b and h — a — 2e lor convenience. Again, the coordinate 
transformations can be obtained, but there is freedom of choice of g and h that can lead 
to more or less convenient coordinate transformations. The linearizability criteria can 
then be given in the form [20] 

3(c/i)^ + Mcy - 2ggx - ghy - Sc^x - '^gxy - hyy = 0, 

S{dg)y + Scdx - 2hhy - hgx - Sdyy - 2hxy - gxx = 0. (9) 

We shall be using this form. 

This equation can be differentiated to yield a third order equation of the general form 

y'" + {A,y'' - A,y' + A,)y" + B,y'' - B,y'' + B^y" - B,y' + i?o = 0, (10) 

which is linear in the second derivative, but with a coefficient quadratic in the first deriva- 
tive, where the coefficients are of this equation are identified with the original equation 

by, 

c = ^2/3,^ = Ai/2,/i = Ao, (11) 
d = - y B2dx + k{y) = J{Bi- Aox)dy + l{x), (12) 
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the constant of integration being fixed by compatibility with the metric tensor, provided 
the compatibihty conditions 

B4 = A2y/3, = - B2 = Aoy - Ai^/2, (13) 

hold and the linearizability conditions (6) are satisfied. Once again, the coordinate trans- 
formations are obtainable and provide the solution. 

The above third order ODE is a total derivative. Though it may not be immediately 
obvious that a given ODE is an exact derivative, the equation in this form may seem 
trivial in some sense. One can use (8) in (10) to replace the second derivative term and 
obtain a third order ODE quintically semi-linear in the first derivative, 

y'" _ ay" + (5y" - ^y'^ + 5y'^ - ey' + = 0, (14) 

where 

a = 3c^, (5 = 5cg + Cy,'y = Ach -\- 2g^ + gy — Cx, S — 3cd + 3gh + hy — gx, (15) 

and the compatibility conditions are 

e = 2dg + h"^ + dy — hx, (f) — dh — dx- (16) 

This equation is not an exact derivative, and hence provides a non-trivial utilization of 
the procedure given. The coefficients can be inverted to obtain the original coefficients 
and hence the metric coefficients, yielding the solution. 

We shall use the above formulation to obtain the linearizability criteria for higher order 
equations in the next section. 



3 Conditional linearizability criteria for fourth order 
ODEs 

On differentiating (14) for the scalar third order ODE, writing the independent variable 
as xand the dependent variable as y, we get 

y"" - {bay" - A/3y'' + 3^y" - 25y' + e)y" - ay'' + i(3y - ax)y" - (7, - f3x)y'' + 

{Sy - lx)y'^ - {ey - 6x)y'^ + {4>x - ex)y' + 0a. = 0. (17) 

The general form of this equation is 

/''_(A4y'4_A3/3+^2/'-Ai|/'+Ao)/-56l/''+i?52/''-i?42/''+532/''-i?2 

(18) 



r. 



subject to the identification of coefficients 

c = /V15, 9^ {As- 4cy)/20c, h={A2- 69^ - ^9y + 3c,) /12c, 

d = (A - - + 25^)/6c, (19) 

witfi tlie constraints 

56 = A4/5, 55 = Asy/A - 54 = A2y/^ - 53 = Aly/2 - 

B2 = Aoy — Aix/2, Bi = dhy + hdy — dxy — Aox/2, Bq = dh^ + hdx — dxx, 

Ao = 29d + h^ + dy-hx. (20) 

In tliis form it is a total derivative and may be regarded as trivial in some sense. 

Replacing the second derivative by the first derivative expressions from (8) we get the 
fourth order ODE semi-linear in the first derivative in the seventh degree 

y"" + Pry" - P,y" + P,y" - PAy" + P,y" - P2y" + Pw' - Po = 0, (21) 

which is not a total derivative. For consistency of the identification of coefficients 

c = (P7/I5)'/', 9 = Pe/SSc" - 2cj,/7c, 
h = P^l21(? - 26^V27c + c^/3c - 9Cyl'i(? - ^9y/21c - Cyy/21c^, 
d = P4/21c^ - 38y/i/21c - 2//7c^ + 8yc^/21c^ - %hcy/2lc^ 

+9x/^c - 99y/^c^ - '2hy/7c + 2c^j,/21c' - ^^/Slc', (22) 

so that all four coefficients are explicitly given. Wc now have four differential constraints 

Pa = 2Scd9 + ISc/i^ + 129'^h - 3{h + d)cx + 4rfcy - (2c/ + 3/i)c/, + (3/i + 2d)9y 
-(c + ^9)hx + 2{9 + h)hy - Scd^ + (c + 29)dy + 9xx - 2hxy + dyy, 
P2 = 18chd + 89^d + 79h^ — 6dcx — 5h9x + 4:d9y — 4:9hx + Ahhy — 3cdx + 29dy 

~^9xx 2hxy ~\~ dyy, 

Pi = 6cd^ — 89hd + — 4:d9x — 3hhx + 3dhy — 29dx + 2hdy + hxx — 2dxy, 

Po = 29d'^ + h^d - 2dhx - hdx + ddy + dxx, (23) 

apart from the earlier stated linearizability conditions (9). 



Instead of starting with the quintically semi-linear third order ODE, we can start with 
the form involving the second derivative (10). In this case we get 

y"" + {A^y'^ - A^y' + Ao)y"' + {B.y' - B,)y"' + {C^y'^ - C^y" + C,y' - Co)y" 

-{D.y"' - D^y'^ + D^y'^ - D^y'^ + D^y' - Do) = 0, (24) 

which is scmilincar involving the third third derivative with a coefficient quadratic in the 
first derivative, second derivative squared with a coefficient linear in the first derivative, 
and otherwise quintic in the first derivative. It is a total derivative with the identification 

c = A2/3 = Pi/6, 9 = A,/2 = So/2, h = Ao, (25) 



the differential constraints 



C3 — 7A2j;/3, C2 — hAiy/2 — 2A2X, Ci — SAoy — 2Aix, 
— A2yy/3, D4 — Alyy/2 — 2A2xy/^, -D3 — ^02/j//2 — Aixy + A2xx/'^, (26) 

and the additional constraints that yield d: 

Cq = dy — 2Aqx, Dq = dxx, (27) 

which gives the value as a single and a double integral to be solved simultaneously 

J{Co + 2Aox)dy + k{x) ^ j{j Dodx)dx + l{y)x + m{y), (28) 

where k,l,m are arbitrary functions of the respective variables. We have the additional 
requirements for Di and D2, 

Di + 2Cox + SAoxx = 0,D2- Coy - Aixx/^ = 0. (29) 

One can now use the third order equation (10) to replace the third derivative term in 
(24) and obtain the equation in a form quadratically semi-linear in the second derivative 
that has a coefficient quadratically semilinear in the first derivative for the higher order, 
quartic in the first derivative for the one linear in the second derivative and otherwise of 
sixth order in the first derivative, 

y"" + iQiy' - Qo)y"' - W - R,y" + R2y'^ - Riy' + Ro)y" 

-iSey"' - S,y"' + S,y'' - Ssy'' + S2y'' - S^y' + ^o) = 0, (30) 

with the identification of the first three coefficients given by 

c = Qi/6, g = go/2, h={R2-Ql + Qix - ^Qoy/'2)/Qi. (31) 

Since the coefficient h is very complicated we shall use it as a shorthand for the above 
expression. The fourth coefficient is given by the solution of 

d^ J{Ro-h'^ + 2hx)dy + k{x). (32) 

In the subsequent equations replacing dy by the integrand above and writing d^ (as there 
is no convenient expression for it), the following constraint equations must hold 

R4 = Q?/4, R3 = QiQo + 7Qiy/6, Ri = 2Qoh + + Shy - 2Qox, 

Se = QlQly/12, S5 = -QlQlx/S6 + QoQly/G + QlQoy/4 + Qlyy/6, 

S4 = -QoQix/6 + hQiy/6 - QiQox/4. + QoQoyf^ + Qihy/2 - Qixy/S + Qoyyf^, 
S3 — Qi{Rq — — hx)/2 — hQi/6 — QqQqx/'^ + hQoy/2 + Qohy + Qixx/^ — Qoxy + hyy, 

S2 — Qidx/2 + Qq{Rq — h'^ + hx) — hhy — hQox/'^ + Roy + Qoxx/'^, 

Si — h{Ro — h'^ + 5hx) — 2Rqx — Qodx — Shxx, Sq — hdx — d^^,(33) 
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along with the hnearization criteria. This is not a total derivative. Though looking more 
messy, it is perfectly usable. 

Yet another form that is not a total derivative can be obtained from (24) by using (8) to 
replace the second derivative term. In this case we get 

y"" + {A^y'^ - A,y' + A,)y"' + B,y" - Bey" + B,y" - B,y" + B,y" 

-B2y'^ + B,y'-Bo^0, (34) 

with the identification of all four coefficients given by 

d= {B4- ScAiAo - ^Al + 3AiCx - TAoCy + 2cAi^ 

-5AiAiy/4: - ScAo^ + 2c^y - ^Aiyy)/4:A2, (35) 

where we have largely used the symbol c rather than A2 to avoid the extra complications 
due to a square root and its differentiation. Further, though d is given algebraically, it is 
too complicated to use conveniently. The constraint equations are 

By = 6c^ Be = 7c{Ai + Cy), B^ = 3^2 + 5cAl/2 - 6cc^ + 7AiCy/2 + 5cAiy/2 + Cyy, 
B3 = ScAid + 6cAl + AjAo - GAqc^ + 7dcy - AiAi^ + 5AoAiy/2 

— 2cAqx + 3AlAQy/2 + Cdy + C^X ~ Al^y + A^yy, 

B2 = 12cAod + Aid + AiAl - Gdc^ - 2AqAxx + hdAxyl2 - A^A^^) 

+3A0A0J, + A\dyl2 + A\xxl2 — 2Aoxy + dyy, 
Bi = 6cd^ + 2AiAod - AdA^^ - 2AqAqx + MAoy + A^dy + Aq^x - 2dxy, 

Bq = d{Aid + dy- Aqx) + 4x,(36) 

along with the linearization conditions. Again, this seems very messy but it has the 
advantage that all four coefficients are identified. 

We could have one other procedure, to use both equations (8) and (10) to replace the 
second and third derivatives. Since this only involves the first derivative to the same 
power as before, this will not yield a new class. Thus we have the following theorems. 

Theorem 1: Equation (21) is linearizable with the identifications (22) if the constraints 
(23) and the linearizability criteria (9) are satisfied, where (c 7^ 0^, with h — a — 2e, g — 
f — 2b, after requiring consistency of the Christoffel symbols with the deduced metric 
coefficients. 

It is to be noted that the four linearizability conditions (6) are stated in terms of the 6 
coefficients a, / and not the four coefficients c, d, g, h. Thus there is degeneracy in the 
choices available. Any choices of a and e for a given h, or /and b for a given g, compatible 
with the metric coefficient relations (7) are permissible. For each such choice we would 



get corresponding linearizability conditions. Instead, we use the conditions (9) to check 
hnearizabihty and exploit the freedom of choice to construct a convenient metric. 

Theorem 2: Equation (30) is linearizable with the identifications (31) and (32) if the 
constraints (33) along with the linearizability criteria (9) are satisfied, after requiring 
consistency of the Christoffel symbols with the deduced metric coefficients. 

Theorem 3: Equation (34) is linearizable with the identifications (35) if the constraints 
(36) and the linearizability criteria (9) are satisfied, after requiring consistency of the 
Christoffel symbols with the deduced metric coefficients. 

Note that there are two other equations, (18) and (24), that are total derivatives of 
linearizable equations with appropriate identifications, that are also linearizable subject 
to the corresponding constraints and conditions. We do not stress on these because they 
may be regarded, in some sense, as trivial. 



4 Examples 

In this section we present some examples of fourth order equations that can be linearized 
by our procedure. 

1. The equation 

is a fourth order total derivative of the form of (18). However, it is clearly not trivial to 
spot this fact by looking at the equation. It satisfies the required constraints (20) and 
linearizability conditions (9). 

2. The equation 

y"» _ 2Wy3 ^ zZky'^ly" + (28i - l^k^)y''' jy + {k^ - 2>hl)ky' /2 - - 5l)ly = 0, (38) 

which is a fourth order equation that is not a total derivative, of the form (21) and is 
linearizable as it satisfies the constraints (23) and the linearizability conditions (9). 

3. The equation 

y"" - {Ay'/y - k)y"' - ^y"^/y + {IQy'^/y' - l)y" - ^y"/y' = 0, (39) 

which is a total derivative of the form of (24) and may be more easily identified as a 
total derivative. It is linearizable as it satisfies the constraints (26) and (29) and the 
linearizability conditions (9). 



n 



4. The equation 



/''-4/Vy'-(6y'Vy'-8V/y+^'+0/+VVz/'-2VVz/'- WVZ/+W = 0, (40) 

which is of the form (30) and satisfies the constraints (33) and the hnearizabihty conditions 
(9). It is not a total derivative. 

5. The equation 

y"" - {^y'y - k)y"' + ^y"/y' + ?,ky''/y^ + (4/ - k')y''/y - Ikly' /2 - Sfy = 0, (41) 

which is not a total derivative and is of the form (34) satisfying the constraints (36) and 
the hnearizabihty criteria (9). As such it is linearizable. 

Though it is not apparent from looking at the equations they come from differentiating 
the linearizable third order equation 

y'" - Gy'^y' + Sky'^/y - {k' - 5l)y' + kly = 0, (42) 

and the total derivative linearizable equation 

y'" - {Ay'/y - k)y" + ^y'^y^ - ly' = 0, (43) 

(in some cases subject to conditions of some equation holding) which both come from 
difi^erentiating the second order linearizable equation 

y"-2y'yy + ky'/2 + ly^0, (44) 

and in one case applying a condition. As such, their solutions necessarily have two arbi- 
trary constants but may not have more. It may be noted that these equations do not have 
any explicit dependence on x. That symmetry being guaranteed only one more needs to 
be looked for to find the solution by symmetry analysis. Using the procedure of writing 
the equation as of geodesic type, we can write down the solution directly. 

The fact that all of these fourth order equations have a common root becomes obvious 
from our analysis right in the beginning, as the identification of the four coefficients in 
each case gives 

c = Q^g = 2,h = k/2,d= -I. (45) 

6. The fourth order ODE 

y"" + {3xy'^ + 2/x)y"' + 6xy'y"^ + {6y'^ - A/x'')y" + Ay'/x^ = 0, (46) 

is of the form of (24) and is hence a total derivative (though this fact is not obvious by 
inspection), with the coefficients given by (25), (27) and (28). It satisfies the constraints 
(26) and (29) and the hnearizabihty criteria (9). As such it is linearizable. 

1 n 



7. The fourth order ODE 



is of the form of (30), and hence is not a total derivative, with the coefficients given by 
(31) and (32) and satisfies the constraints (33) and the hnearizabihty criteria (9). As such 
it is hnearizable. 

8. The fourth order ODE 

y"" + (3xy'^ + 2/x)y'" + 6x^y'^ + 18y'^ + 16y'Vx + Vly' jx^ = 0, (48) 

is of the form of (34), and hence is not a total derivative, with the coefficients given by 
(35) and satisfying the constraints (36) and the hnearizabihty criteria (9). As such it is 
hnearizable. 

9. The fourth order ODE 

y"" - {Ihx^y'^ + 21?/'2 + 6/x')?/" - 6xy'^ + Vly' jx"- = 0, (49) 

is of the form of (18), and hence is a total derivative, with the coefficients given by (19), 
satisfying the constraints (20) and the hnearizabihty criteria (9). As such it is hnearizable. 

10. The fourth order ODE 

y"" + (3x|/'2 + 2lx)y"' + ^x^y" + \%y'^ + l^y'^-'jx + 12y7a;^ = 0, (50) 

is of the form of (21), and hence is not a total derivative, with the coefficients given by 
(22) and satisfying the constraints (23) and the hnearizabihty criteria (9). As such it is 
hnearizable. 

It is again apparent that even the total derivative equation is not obviously so, and that 
the other equations being hnearizable would not be clear by inspection. The coefficients 
of examples (6) to (10) are all c = a;, /i = 2/x^g = c? = and the root equation is the 
geodesic equation for flat space in polar coordinates. Thus the linearizing transformation 
is simply the conversion from polar to Cartesian coordinates. 

It will be noticed that in the above examples either x or y are missing from the coefficients. 
This makes the application of the identification and constraints relatively trivial. We end 
with a couple of non-trivial examples, in which both arise. 

11. The fourth order ODE 

y"" - 15xy^/y^ - Ibxy'^y^ + 397/'Vl/^ + W^/xy^ - ^Qy'Vxy^ - 36y'^/x^y + 24:y'/x^ = 0, 

(51) 

is of the form of (21), and hence is not a total derivative. Its coefficients are given by 
(22) and come out to be c = —x/y'^, g — l/y,h — 2/x,d — 0. They satisfy the constraints 

1 1 



(23) and linearizability criteria (9). As such this equation is hnearizable. In fact it can 
be derived by differentiating the hnearizable third order equation [19] 



(52) 



writing the root equation from the coefficients, and using it to substitute the y" arising 
from the differentiation. The solution is 



y"" - {Qxy'/y' + 2/y)y"^ - {S^x^^jy^ - 2xy''/y' - 7y'^/y^ - Sy'/xy + 8/x^)/ 

+Qxy'/y'> - 2xy">/y' - 2y"/y' - Qy'^/xy^ - Ay'^/x^y + Sy'/x' = 0, (54) 



is of the form of (30). Its coefficients are given by (31) and (32) and are the same as of 
the above example. Hence they share a common root equation. This equation arises by 
differentiating the root equation twice and using its first derivative to replace the y'" . It 
satisfies the constraints (33) and the linearizability criteria (9) and is hnearizable, yielding 
the same linear equation and possessing the same solution as of the previous example. 



5 Concluding remarks 

We have provided a procedure to determine the linearizability of some fourth order scalar 
semi-linear odes that are hnearizable. We have written them as five classes. There could 
have been other ways of getting to the five classes, e.g. by first differentiating and then 
replacing or first replacing and then differentiating. Since, the general form of both 
procedures would be the same they are not different. Again, we could have used a lower 
order equation to partially replace terms in the higher order equation. For example, we 
could have used the second order equation to replace the quadratic term in the second 
derivative in (52) but retain the linear term as it is. These are not really independent 
in some sense. In this sense there are five independent classes. Two of these may be 
regarded as trivial as they are total derivatives. However, a glance at the equations 
will show that even these are not so easy to identify as total derivatives. The other 
three classes are not total derivatives. All can be thought of as arising from the same 
second order hnearizable differential equation that satisfies the Lie conditions, by double 
differentiation, but conditional to the original differential equation and in some cases the 
differentiated differential equation. As such, the linearizability classes are non- classical 
and would not lie in the three classical classes reducible to the fourth order odes y"{x) = 0, 
y"{x) + y{x) = and y" + a{x)y{x) = 0. These are guaranteed four arbitrary constants 
while our solutions are only guaranteed two. 



Axy + Bx/y = 1, 



(53) 



where A and B are arbitrary constant real numbers. 
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Let us call the underlying second order equation the root equation, and the specific fourth 
order equations forms of the similar fourth order equations. It is clear that all forms 
arising from a common second order equation will form an equivalence class. As such, 
all similar fourth order equations will have two solutions in common but may have other, 
different solutions. The question arises whether these equivalence classes of linearizable 
fourth order semi- linear odes are disjoint and can be used to decompose the space of those 
odes that are linearizable of this type. In this context, it seems unlikely that this be the 
case because there would be some of them that are linearizable and have two solutions 
in common with others. As such, we conjecture that the set of linearizable fourth order 
odes with common root equations will not be decomposable into disjoint classes. 

It is worth noting that there is only one class of linearizable second order semi-linear ode 
(as proved by Lie) and two classes of third order odes, one of which is a total derivative. 
For the fourth order we find five of which two are total derivatives and three are not. Can 
this procedure be carried further? The answer is "obviously it can" . The fifth order will 
have one form that only involves the first derivatives apart from the fifth derivative and is 
of ninth order in them. It is obvious that here all four coefficients will be easily identified, 
but now instead of the four constraints (22) we will have six constraints. Similarly, for 
the sixth order the corresponding class will have the first derivative to the eleventh power 
and will have eight constraints and so on. Of course the usual linearizability conditions 
(6) would also have to be met. 

The question arises as to the number of classes that wc can now have. Clearly, for the 
fifth order there will be five classes of total derivative forms obtained by differentiating 
the five classes of the fourth order. How many other classes will there be? The number 
can be obtained by counting the different forms that one could obtain and turns out to 
be six. Thus the total number of classes is eleven. It is obvious that there will then be 
eleven classes of total derivative sixth order equations obtained. One could continue the 
procedure of counting but it would be nice to have a general formula for the total number 
of classes for any order. 

It is worth stressing that as can be seen from the examples, our procedure provides the 
solutions of the linearizable semi-linear equations of our classes for any order that is 
identified. It should not be too difficult to prepare an algebraic code for identification 
of coefficients, the constraints and the linearizability conditions for the various classes. 
Certainly, the form involving only the first order could be relatively trivially provided. It 
would be worth while to get the algebraic codes prepared for at least some of the higher 
order classes. 
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